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Excitation of low-frequency QPOs in black hole accretion 
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Abstract. We study possible mechanisms of excitation of quasiperiodic oscillations in the accretion flow of black 
hole accreters in their hard spectral states, in the context of the 'truncated disk' model. Quasi-spherical oscillations 
of the inner ion-supported accretion flow (ISAF) can be excited by the interaction of this hot flow with the cool 
disk extending outward from it. The fundamental mode of (p-mode) oscillation is most easily excited, and has 
a frequency near the Kepler frequency at the inner edge of the cool disk. The strongest excitation mechanism is 
a feedback loop involving cooling of the ISAF by soft photons from the cool disk and heating of the cool disk 
by the ISAF, while synchrotron emission can be a relatively strong damping effect. Growth times are computed 
by detailed Comptonization calculations of the interaction of an idealized ISAF with a cool disk. Typical growth 
times as short as a few dynamical times are found, while amplitudes can be as large as 10%. 
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1. Introduction 

A common feature of the X-ray light curves of black 
hole binaries (most of them transient X-ray sources) is 
quasiperiodic oscillations (QPOs) (narrow power peaks 
with high quality factor Q > 2 (i.e. Belloni et al. 
2002)). QPOs can be categorized -according to their peak 
frequencies- as high-frequency QPOs (HFQPOs) and low- 
frequency QPOs (LFQPOs). HFQPOs appear at i/>60 Hz 
and have low rms amplitudes (typically 1% of the mean 
count rate). They are typically observed when black hole 
binaries are in the very high state, i.e. when they exhibit 
both a thermal component at fc^T ~ 1 keV and a hard X- 
ray power law with photon energy index (dN/dE oc E~ r ) 
r ~ 2.5. Because of their high frequencies, they are 
thought to convey information about physical properties 
that occur near the black hole horizon. A model which 
explains the 2:3 frequency ratio often seen in HFQPOs is 
given in Rezzolla et al. (2003). 

The more common LFQPOs have a characteristic fre- 
quency range of 0.1-10 Hz, contribute significantly to 
the total variability (rms amplitudes of the order 10% 
and higher have been seen in XTE J1550-564 and XTE 
J1118+480) and have been observed in both the low/hard 
and very high spectral states (for a review see Swank 
2001). Their amplitude is generally higher at higher pho- 
ton energies, though exceptions occur (for example in a 
QPO seen in XTE J1550-564 (Kalemci et al. 2001) and in 
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Cygnus X-l (Rutledge et al. 1999)). Characteristic for the 
low frequency QPOs is a decrease of the frequency during 
the decay of the transient, suggesting a physical relation 
between frequency and accretion rate. 

In recent years, it has become increasingly clear that 
the hard state power density spectrum (PDS) of Cygnus 
X-l and other black hole binaries shows several breaks, i.e. 
frequencies where the PDS steepens (Nowak et al. 1999) 
and that this complex structure can be described with 
multiple broad noise components (Nowak 2000; Nowak 
et al. 2002; Pottschmidt et al. 2003). Moreover, there is 
evidence of a correlation between the horizontal branch 
oscillations and lower kHz QPOs in Z and atoll sources 
that can be extrapolated, by more than two orders of 
magnitude in frequency, to include a correlation between 
LFQPOs and noise components in black hole binaries 
(Psaltis et al. 1999). In this work, we develop a model 
that attempts to clarify the origin of at least some of the 
LFQPOs. 

1.1. QPOs as diagnostics of the accretion flow 

Apart from their intrinsic interest, the oscillations are a 
potentially important diagnostic of the still poorly un- 
derstood inner regions of the accretion flow. In the 'disk 
corona' model (e.g. Poutanen & Fabian 1999), a cool (< 1 
keV) disk extends to the innermost stable orbit around the 
hole, while the hard X-rays are produced in a magnetically 
heated corona above it. In the 'truncated disk' model, on 



2 



Dimitrios Giannios and Henk C. Spruit: Excitation of low- frequency QPOs 



the other hand, the cool disk extends inward only to a 
transition radius R tI , inside of which the accretion takes 
place as an ion supported flow (Rees et al. 1982). In such 
flows the ion temperature is near the virial temperature, 
while the electrons producing the hard X-rays have a tem- 
perature around 100 keV (Shapiro et al. 1976; Ichimaru 
1977). Their hydrodynamical properties have been stud- 
ied by Gilham (1981) and in detail by Narayan and Yi 
(1994 and references therein) who called them ADAFs. 
To distinguish them from the radiatively supported accre- 
tion flows (also called ADAF), we use here the name 'ion 
supported accretion flow' (ISAF). 

The arguments for and against either of these geome- 
tries are part of a current debate in the literature. The 
main obstacle to acceptance of the truncated disk model 
is the difficulty in visualizing a process that can make a 
cool disk flow to change into a hot ion supported flow, 
across a highly unstable region of intermediate tempera- 
tures. Without going into details, we note here that this 
difficulty has been overcome by Spruit & Deufel (2002). It 
is shown there that viscous heating of the ions in the in- 
teraction region between the ISAF and the cool disk near 
the transition radius causes the cool disk to evaporate and 
feed the ISAF even when the transition radius is close to 
the last stable orbit. 

Distinguishing the two accretion geometries observa- 
tionally is not possible from the X-ray spectrum alone. 
The hope is that a prominent feature like the low- 
frequency QPO can help in settling the question of the 
accretion geometry. In the truncated disk model, the trun- 
cation radius provides a characteristic length scale and 
an (orbital) frequency scale that can be associated with 
the QPO. The dependence of frequency on accretion rate 
would then imply an increase of the transition radius with 
decreasing accretion rate. Supporting evidence for such a 
physical relation is found in the correlation between re- 
flection indicators (Fe line strength and amplitude of the 
Compton reflection hump) and the strength of the soft 
component in the X-ray spectrum with QPO frequency 
(Zdziarski et al. 1999; Gilfanov et al. 2000). 

1.2. The importance of excitation 

The number of possible modes of oscillation in an accre- 
tion disk is as large as those of a star. Discrete spectra of 
trapped oscillations can occur near the last stable orbit 
around a hole because of the shape of the relativistic po- 
tential (Kato 1990; Honma et al. 1992; Ipser 1996; Nowak 
& Wagoner 1991, 1992). But even a disk in a Newtonian 
potential supports trapped oscillations near its inner or 
outer edge (Drury 1980; Blaes & Hawley 1988; Kaisig 
1989). The spectra of such modes are quite complex (e.g. 
Narayan et al. 1987). It is not surprising, therefore, that 
picking modes to match an observation is not too difficult 
(e.g. Titarchuk et al. 1998). 

To make progress in interpreting QPOs, criteria are 
needed to identify the most plausible modes. An impor- 



tant clue is the amplitude of the oscillation. To explain 
modulations in the total X-ray luminosity as high as 10% 
with a hydrodynamic oscillation, the amplitude of the os- 
cillation must be significant, not only at the larger dis- 
tances where the frequency of the oscillation is presumably 
set, but also in the regions close to the hole where most 
of the X-rays are produced. Secondly, the observed flux is 
an average over the accretion flow, so that all nonaxisym- 
metric modes average out (except for strong nonlinearity, 
relativistic beaming, or optical depth effects). Similarly, 
the visibility of modes with one or more nodes in the ra- 
dial part of their structure will be strongly suppressed. 
The most likely modes capable of producing large ampli- 
tudes in the X-ray flux are therefore axisymmetric modes 
with few nodes in the radial and vertical directions. 

Finally, there must be a plausible mechanism to excite 
the mode. While broad band noise might be explainable 
as some consequence of turbulence in the flow (Lyubarskii 
1997; Churazov et al. 2001), the high amplitudes seen in 
LFQPOs probably require an excitation process intrinsic 
to the oscillation itself. 

The wide and possibly continuous variation of the 
LFQPO with (inferred) accretion rate is natural in the 
truncated disk model. In addition, the presence of two 
components in the accretion flow around the transition 
radius gives rise to possible excitation mechanisms that 
are absent in disk+corona models. 

2. Excitation by interaction 

An ion supported accretion flow and a cool disk partially 
embedded in it interact through exchange of mass, energy 
and angular momentum. Exchange of mass takes place in 
the form of an evaporation process from some part of the 
cool disk, as needed to produce the ISAF in the first place, 
and possible condensation onto the disk at other locations. 
The evaporation can take place at a large distance where 
the virial temperature is of the order 10 6 -10 8 K, through 
processes similar to those maintaining the mass balance in 
the solar corona (Meyer et al. 2000; Meyer-Hofmeister & 
Meyer 2001; Liu et al. 2002). The observations mentioned 
above, however, show that the cool disk is truncated at a 
much closer distance to the hole, at least in some cases. 
The mechanism that allows the disk to evaporate at such 
distances has been identified by Spruit & Deufel (2002). 

Thus assuming that the cool disk evaporates from a 
region near its inner edge, the resulting ISAF spreads 
both inward, feeding the hole and outward, spreading over 
the cool disk. At some distance from the inner edge, the 
spreading ISAF may condense back onto the cool disk, so 
that a region of coexistence between the cool disk and the 
ISAF is maintained. 

An oscillation might be excited by the exchange of 
mass between the cool disk and the ISAF because the 
accretion time scale in the disk is much longer than the 
viscous spreading time of the ISAF. Mass evaporating 
from the truncation radius R tr depletes the surface den- 
sity there. When this mass condenses from the ISAF again 
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on to the disk at some distance i? ut > -Rtr, it cre- 
ates an enhancement there which takes a viscous time 
(R ou t — Rtt) 2 /v to arrive at R tr - Spreading in the ISAF 
and condensation at i? ou t can then feed the next circle. 

A second possibility for excitation is through inter- 
action by angular momentum. A known example is vis- 
cously excited oscillation found by Kato (1978), and stud- 
ied in detail by Honma et al. (1992); Papaloizou & Stanley 
(1986) and others. In this mechanism, a sound wave oscil- 
lation is excited viscously in a rotating flow. The ampli- 
tude of the oscillation peaks at the outer edge of the flow 
(Kley et al. 1993). The boundary between the ISAF and 
the cool disk might be sufficient to excite a sound wave in 
the ISAF through this mechanism. 

Finally, exchange of energy between the ISAF and the 
cool disk can excite oscillations. The flow is cooled by 
soft photons (cither internally generated by synchrotron 
emission or externally from the cool disk) and heated by 
viscous dissipation. Changes in the soft photon input heat 
or cool the ISAF, causing it to expand and contract on 
a sound crossing time scale. Since the ISAF is partially 
pressure supported, the relative changes in size are of the 
same order as the relative changes in heating or cooling 
rate. The size of the region of overlap with the cool disk 
thus varies as well, and thereby the soft photon input into 
the ISAF. 

In the following we investigate the conditions under 
which interaction by exchange of energy can produce self- 
excited oscillations. More specifically, we focus on the vari- 
ation of radiative cooling of the flow, and assume that the 
viscous heating does not vary. The viscous heating can ex- 
cite/damp oscillations on the viscous time scale. However, 
we show in section 5 that radiative cooling excites/damps 
oscillations on much shorter time scales, validating our 
assumption (for further discussion see Sect. 5.3). 

3. Model 

We consider low order oscillation modes of the ISAF, ex- 
cited by the energetic interaction between the ISAF and 
the cool disk. There are different interaction paths, and 
whether these excite or damp an oscillation of the ISAF 
depends on geometrical properties and on the source of 
soft photons. 



a black hole if the magnetic field strength is some frac- 
tion of the equipartition field. This synchrotron radiation 
is strongly self absorbed below a characteristic turn-over 
frequency v t ■ Most of the synchrotron flux is emitted near 
this frequency, which typically lies in the infra-red or the 
optical region of the electromagnetic spectrum. 

A second source of soft photons is the radiative in- 
teraction between the hot flow and the cool disk. Hard 
X-rays from the inner hot flow are partly absorbed and 
partly reflected by the cool disk. The absorbed energy is 
reprocessed into thermal soft flux and emitted by the disk. 
These soft photons are the seed photons Comptonizcd 
by the ISAF. The energetics and spectra produced in 
this radiation-mediated interaction have been studied by 
Haardt & Maraschi (1993). 

In addition to this purely radiative interaction, the 
hot flow also heats the cool disk by 'proton illumination' 
(Deufel et al. 2001; Deufel et al. 2002; Spruit & Deufel 
2002). In this process, part of the ISAF condenses onto 
the disk in the overlap region. Energetic protons (which 
carry most of the ISAF's energy) are slowed down by the 
disk plasma and their energy gets thermalized and appears 
as soft flux. The proton energy flux F p which is released 
in the disk is of the order 



F p = n p (k B T p )c s 



(1) 



where c s is the proton sound speed and n p is the proton 
number density in the ISAF. The density is of the order 
n p = t/(<jt Rout), where r is the Thomson optical depth 
of the ISAF and <tt the Thomson cross section. 

If half of this flux is re-emitted isotropically as soft 
photons, the disk surface temperature Td; s k is given by 
the relation F p /2 = crT^ isk , with a the Stefan-Boltzmann 
constant (a — 5.67 • 10~ 5 erg cm~ 2 s _1 K~ 4 ). We get 

T disk * (^7j) 1/4 keV, (2) 

where mi is the mass of the black hole in units of 10 solar 
masses (2 • 10 34 g) and r the distance from the black hole in 
units of the Schwarzschild radius R s = 2GM/c 2 . We see 
that 7di s k ~ 100 eV is a typical value for r of order unity 
and r ~ 50. In the calculations reported below, the effects 
of each of these sources of soft photons on the oscillation 
are investigated. 



3.1. Geometry and soft photon sources 

The basic ingredients of the model are a standard thin 
disk (Shakura & Sunyaev 1973) truncated at a radius R tI 
and an ISAF extending to a radius i? ut > Rtv Inverse 
Compton scattering of soft photons is the natural radiative 
process to explain the non-thermal X-ray spectrum in the 
low/hard state, but there are several possible sources for 
the soft photons, depending on the geometry of the flow 
and model assumptions. 

Soft photons can be generated internally in the ISAF 
as thermal synchrotron photons. This will be the domi- 
nant source of soft photons in an ion supported flow near 



3.2. The oscillations, model approximations 

The excitation mechanisms considered here operate on the 
balance between heating and cooling in the ISAF, in par- 
ticular through a modulation of the soft photon input. 
The simplest modes and the most likely ones to be ex- 
cited by this mechanism are the axisymmctric modes with 
few nodes in the radial direction, in particular the funda- 
mental mode. In the following, we restrict ourselves to the 
fundamental mode of the ISAF. As in a star, it has a fre- 
quency of the order of the Kepler frequency at its outer 
edge, or the inverse time for a sound wave to travel across 
it. 
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Since the fundamental mode senses only the average 
sound speed, it is sufficient to model the ISAF as spherical. 
Since the rotation period is of the order of the sound cross- 
ing time, its effects on the fundamental mode are modest, 
hence we also ignore the effect of rotation on the mode. 
Each of these approximations can be easily improved on 
with a more detailed model of the structure of an ISAF. 
At the present level of understanding of the geometry of 
the accretion flow onto black holes, this is not a primary 
concern. 

We further simplify the model by taking the equilib- 
rium electron temperature T c and density po as homoge- 
neous. The fundamental mode of oscillation is approxi- 
mated by a homologous expansion and contraction. 

Thus, if i?(t, Rq) is the radial position of a shell of mass 
with equilibrium position Rq, a sinusoidal oscillation of the 
model ISAF is described by 



R(t, R ) = i? [l + a • sin (uit)), 



(3) 



where a is the fractional amplitude of the oscillation. The 
radial velocity and the density perturbation are 

v(t, i?o) = aujRo cos(u>t), Ap = — 3apo sin(wt). (4) 

Also the perturbation of the optical depth of the ISAF 

(r = p K es Rout) is 



At = —2aTsm(ujt). 



(5) 



Such a non-rotating, constant density configuration in 
homologous motion has a fundamental frequency of oscil- 
lation (Cox 1980) 



w = a/(37-4) 



GM 

^out 



(6) 



where 7 is the ratio of specific heats, assumed constant. 
For the ion cloud of an ISAF, the non-relativistic value 
7 = 5/3 is a good approximation, which will be used in 
this work. The oscillation frequency is then exactly equal 
to the Kepler frequency ujk ■ (In view of the approximations 
made, this exact agreement is a coincidence.) 

3.3. Excitation and damping 

An efficient excitation mechanism is needed, since the ob- 
served oscillation is often quite strong, accounting for a 
large fraction of rms variability. In the present model, the 
excitation is due imbalance of cooling and heating terms 
during the oscillation, such that the net work done on the 
gas by the pressure forces is positive during an oscillation 
cycle. 

With standard stellar oscillation theory (see for exam- 
ple Cox 1980; Unno et al. 1979), the rate of excitation of 
an oscillation can be expressed by the work integral 

f-T 



(W) = - 



dm 



(7- l)^Aqdt. 
P 



(7) 



Where Aq is the (Lagrangian) variation of the net cooling 
rate per unit mass during the oscillation, the first integral 



is over the mass of the ISAF and the second over the 
oscillatory period T. An ideal gas with constant ratio of 
specific heats 7 has again been assumed. 

If the excitation rate is small compared with the oscil- 
lation frequency (as needed in order for the quality factor 
Q of the oscillation to be high), the work function can 
be represented by difference in the net cooling rate at the 
maxima and minima of the oscillation (see Appendix A) 



(W) = 3(7 - l)et 



Qcxp ycontr 



(8) 



where a is the relative displacement amplitude as before, 
and <?oxp, 9contr are the net cooling rates at maximum ex- 
pansion and maximum contraction, respectively. 

If cooling mechanisms are more efficient in the ex- 
panded phase of the oscillation than in the contracted 
phase, then in every oscillatory cycle there is a net increase 
of the oscillation's energy proportional to g oxp — (fcontr 
which means excitation of the oscillation. If, on the other 
hand, q cxp — q con t Y < then the oscillation gradually de- 
cays. 

3.4. Effect of different sources of soft photons 

To find out whether a mode can be excited by the pro- 
cesses that determine the energy balance, special attention 
must be paid to the cooling mechanisms of the ISAF, since 
these are more easily modified than the viscous heating of 
the flow (cf. Sect. 5.3). Inverse Compton accounts for most 
of the cooling of the ISAF (at least for the conditions un- 
der which the observed spectra are produced, since these 
are dominated by hard X-rays). Soft flux, that crosses the 
hot and optically thin medium of the ISAF, is upscat- 
tered by energetic electrons resulting in energy transfer 
from the electrons to the radiation field. The cooling rate 
of the ISAF is thus controlled by the input of soft pho- 
tons. As discussed above, we consider three sources of soft 
photons: 

i) synchrotron photons from the inner parts of the flow, 

ii) irradiation of the thin disk by X-rays from the ISAF, 

iii) proton illumination of the thin disk in the overlap 
area. 

The strength of Compton cooling is determined by 
the parameter, y — Ak^T c / '(m c 2 )r and the soft flux: 
q c = L so f t (e y — 1). During the expansion phase the op- 
tical depth of the ISAF drops (Eq. 0) and so does the 
electron temperature (see next section), resulting in a de- 
crease of the y parameter. The opposite trend is observed 
during contraction. So, if the soft flux does not vary much, 
the cooling rate of the expanded phase must always be less 
than the cooling rate of the contracted phase. This leads 
to damping of the oscillations. On the other hand if the 
soft flux increases quickly enough during expansion, there 
can be excitation of the oscillations. 

When synchrotron photons are the 'seed' photons to 
be Comptonized, then it is straightforward to show that 
they act as a damping mechanism. This is because the syn- 
chrotron flux increases strongly with electron temperature 
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and the strength of the magnetic field and less strongly 
with the Thompson optical depth (e.g. Wardzihski & 
Zdziarski 2000). The magnetic field scales as B 2 ~ P gas , 
which means that B, T c and r all decrease during expan- 
sion. The synchrotron soft flux decreases and so does the 
cooling rate during expansion, resulting in energy loss of 
the oscillation. 

On the other hand, in the irradiation and proton il- 
lumination cases, it is not a priori clear whether the net 
effect is damping or exciting. Detailed calculations of the 
cooling rates are needed to settle this question. 

4. Calculations 

4.1. The electron temperature in equilibrium 

The energy source of the electrons in the two tempera- 
ture plasma is the transfer of energy from the hot protons 
through Coulomb interactions. This causes the protons to 
cool at a time scale (Spitzer 1956; Stepney 1983) t p 



1 



! l/ 2 m c nebcin A \ m c c 



/fc B T c \3/2 

Vm r c 2 / ' 



(9) 



where T c stands for the electron temperature, n = 
t I (ox-Rout) is the number density of electrons, In A is the 
Coulomb logarithm. This expression is the non-relativistic 
one and applies when T <C T p <C m p T c /m c . Where T p is 
the proton temperature assumed to be the virial temper- 
ature 



_ GMm p 156 
p ™~ 3fc B i? " r 



MeV. 



(10) 



As will be shown in the rest of the section, these inequal- 
ities hold for most of the parameter space of T c , T p in 
which we are interested in this work. 

Because of their lower temperature, the heating time 
scale of the electrons is shorter than the proton cooling 
time scale 



_ T e 
t c — 7^rt p . 



(11) 



The heating rate of the electrons (erg/sec) over the whole 
ISAF is then given by nk B T c V/t c , with V being the vol- 
ume of the ISAF. In equilibrium this quantity equals the 
cooling rate of the electrons, which is assumed to be con- 
trolled by Compton cooling and ultimately equal to the 
observed X-ray flux of the source 



nk B T c V 



(12) 



Combining this equation with Eq. (JTTJ, Eq. ©, V = 

(4/3)7ri?Q Ut and taking In A ~ 20, T p — T vil , we can relate 
the temperature to the optical depth and X-ray flux 



k B T c ~ 5r 



4/3 



Le 



2/3 



keV 



(13) 
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Fig. 1. Electron temperature in an ion supported flow 
of optical depth r, as a function of the observed X-ray 
luminosity L x and luminosity index s of the spectrum 
(yF v ~v s ). 



A further constraint on T e and r from observable quan- 
tities comes from the slope s of the X-ray spectrum in 
a v ■ F v diagram (y ■ F v ~ v s ). With an analytic fit 
to Comptonized spectra given in Wardzihski & Zdziarski 
(2000) their Eq. [22], which agrees with reasonable accu- 
racy with our computational results in Sect. 15. 4[) we get 
for the hard X-ray power-law slope 



0.8 lOOkeV 

r 4 / 5 k B T r , 



, 0.5 < t < 2, k B T c < 200 keV. (14) 



(Shapiro et al. 1976). Here Le is the Eddington luminosity. 



Eqs. i|13[l an d (|14|) determine the electron temperature 
of the ISAF, expected from electron heating through 
Coulomb interaction with virial protons, in terms of the 
optical depth t of the flow and two observables, the X-ray 
luminosity and the slope of the hard X-ray spectrum. This 
is shown in Fig. 1. 

From Fig. 1 one can see that the electron temperature 
of the ISAF, if it is to account for the observations, is 
^100 keV. Furthermore since T p = T VU: the inequalities 
T c < T p < m p T c /m c hold for r out < 1000 validating the 
use of Eq. © ■ 

If the optical depth is expressed in terms of the accre- 
tion rate, with an assumed viscosity, this can be further 
reduced to a dependence of the expected electron temper- 
ature on accretion rate (e.g. Rees et al. 1982; Narayan & 
Yi 1995; Zdziarski 1998). 

4.2. The electron temperature during oscillation 

In order to calculate the excitation of an oscillation, the 
net cooling/heating of the flow has to be followed during 
a cycle (Sect. 3.3). For this we need the variation of the 
electron temperature during the oscillation. We compare 
the electron heating time scale with half of the oscillation 
period P (time between two passes through the equilib- 
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rium position). From equations ©, (|10fl . (|llfl the heating 
time scale of the electrons is 



7 mir2 / k B T c \5/2 



sec, 



and 



P/2 = 4.3- 10" 4 m! 



-& 3/2 



sec, 



(15) 
(16) 



where r ou t is the radius of the ISAF in R s . It follows that 
t c < P/2 for optical depths > 0.01 and r out < 1000. 

The electron temperature therefore adjusts quasistat- 
ically to the changing proton temperature, so equilibrium 
equations from the last section still apply during the oscil- 
lation. Taking the Lagrangian derivative (A) of Eq. (|ll|l 
and using An — —3ansm(uit) and AT p = — 2T p sm(uit) 
(see next section) we find that 



5t c [a sin(ujt) + 



2T n 



(17) 



Since the adjustment time of the electron temperature is 
short, the heating rate of the electrons must equal the 
(Compton) cooling rate g C ooi at any instant: 



_ nk B T c V _ . 

9hcat T 9cool ■ 



4.3. Evaluation of the work integral 



(18) 



The oscillations involve the pressure force and gravity. In 
an ion supported flow, both these forces are concentrated 
in the protons. Hence the cooling/heating rate q in the 
work integral involves only the proton temperature and 
the distance from the central mass 



q = q(T p ,r(t)). 



(19) 



The proton temperature T p of a mass element at distance 
r(t) changes with time due to three factors, the compres- 
sion/expansion of the gas, energy transfer to the electrons, 
and the viscous heating. Variations in viscous heating as 
a possible excitation mechanism of oscillations in accre- 
tion disks has been studied before (Kato 1978; Papaloizou 
& Stanley 1986; Kley et al. 1993). The maximum growth 
rate due to this process is of order af2, where a is the disk 
viscosity parameter. Assuming a to be small, perhaps of 
order 0.01-0.1 as indicated by numerical simulations of 
magnetic accretion disk turbulence (e.g. Brandenburg et 
al. 1997; Armitage 1998), this contribution is small com- 
pared with the growth rate we find below due to the other 
interaction processes, hence we will neglect variations in 
viscous dissipation during the oscillation cycle in the fol- 
lowing (see also Sect. 15. 3|) . 

Assuming we are looking at a slowly damped or grow- 
ing mode, the temperature variation can be written as 
AT = AT ad + eAT cool , where AT ad (i) is the adiabatic 
temperature variation, AT coo i is the temperature change 
due to the variation in energy transferred to the electrons, 
and e a small number. Expanding in e, the work integral 
can be evaluated by setting AT = AT a d in Eq. (fP9"|) . that 



is, it can be found to lowest order using only the adiabatic 
variation of the proton temperature. 

In practice we evaluate the work integral in a highly 
discretized way: the integral is represented by the sum 
of two values: the integrands at the phases of maximum 
contraction and maximum expansion. For a linear (sinu- 
soidal) oscillation, evaluation at a single known phase in 
the oscillation would be enough to determine its ampli- 
tude, and would be sufficient to determine the work inte- 
gral. By evaluation at the phases of both maximum and 
minimum expansion, nonlincarity due to asymmetry be- 
tween the expanded and contracted phases can be taken 
into account to lowest order. This is somewhat important 
for the present problem because an asymmetry is in fact 
expected since, for large oscillatory amplitude, the ISAF 
in its contracted phase can actually detach from the cool 
disk. For the two phases, the cooling rate is evaluated by 
an iterative procedure, varying the electron temperature 
until the balance Eq. (|18Jl is satisfied. 

4.4. The Monte-Carlo code 

In order to study both the excitation and the damping 
mechanisms that take place in the ISAF quantitatively, 
the Comptonization of soft radiation has to be calculated 
in detail. For this we have performed Monte Carlo sim- 
ulations following Pozdnyakov et al. (1983). For specific 
values of the electron temperature T e , the optical depth 
r, location and energy of the soft photons, we calculate 
the cooling rate of the ISAF as well as the X-ray spec- 
trum produced. 

The location of the soft photons depends on the mech- 
anism that produces them. In the case of thermal syn- 
chrotron photons, they are expected to be emitted in the 
region close to the black hole, with an energy distribution 
that strongly peaks at a few eV. For simplicity we assume 
that all synchrotron photons come from the center of the 
ISAF. If the soft photon flux is a result of disk irradi- 
ation caused by hard X-rays coming from the inner hot 
flow, then the location of their emission and their energy 
is calculated by the iterative method that is described in 
section 5.1. Finally in the 'proton illumination' case, the 
soft flux comes from the inner part of the thin disk where 
the disk interacts with the ISAF (at radii i?tr < r < -Rout)- 
The photon energy in this case is given by Eq. J2J. 

The Monte Carlo method follows the path of a large 
number of photons and makes a statistical description of 
the radiation transfer inside the scattering medium. The 
calculation takes properly into account the quantum me- 
chanical and special relativistic effects of the microphysics 
of the scattering. We make sure that the statistical sam- 
ple (i.e. the number of photons) is large enough to keep 
statistical noise at a low level. 

5. Results 

Whether small amplitude radial oscillations of the ISAF, 
once appearing, will increase in amplitude or die out de- 
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pends on the origin of the soft photons. As discussed in 
section lXTl there are three possible sources of soft photons: 
the reprocessing of hard (inverse Compton) photons in- 
tercepted by the cool disk (the 'X-ray irradiation' model), 
photons produced by ion-heating of the cool disk ('proton 
illumination'), and internal (synchrotron) photons gener- 
ated in the ISAF. The conditions for excitation are differ- 
ent in each case. Although all these processes may con- 
tribute, we study their effects separately. We start with 
the irradiation model. 

5.1. X-ray irradiated cool disk 

Parameters of the ISAF are its equilibrium values of elec- 
tron temperature and optical depth, using Fig. 1 as a 
guide. The degree of overlap of the ISAF over the cool 
disk in the equilibrium state, / = R ut/Rtr — 1, is a third 
parameter of the problem. We make a guess of the initial 
radial temperature profile of the disk surface and simu- 
late the emission of cool photons from the disk and the 
upscattering of these soft photons into the ISAF with a 
Monte Carlo calculation. Part of the hard X-rays produced 
by inverse Compton scattering in the ISAF will meet the 
thin disk in their scattering random walk and get either 
absorbed or reflected. The absorbed X-ray flux, assumed 
to be thermalized into a black body, determines the new 
temperature profile of the disk surface. The simulation is 
repeated using the new temperature profile several times 
until it converges (within numerical uncertainties imposed 
by the statistical nature of the Monte Carlo method) . The 
cooling rate of the ISAF in the final run is recorded. 

The next step is to expand the medium by a small 
amount. The optical depth of the hot medium is reduced 
in accordance to Eq. J^J for sin(wt) = 1 (maximum ex- 
pansion), while the temperature of the hot medium on the 
expanded phase is given a trial value. Following the proce- 
dure described in the previous paragraph, the temperature 
profile of the disk and the cooling rate of the ISAF q exp 
are computed. Then we iterate varying the temperature of 
the ISAF until Eq. Ijl8(l is satisfied. The same procedure 
can be applied to maximum contraction of the ISAF to 
compute its cooling rate g con tr- 

There are two simplifications made during this com- 
putation. First, we have ignored the reflected component. 
Detailed spectral predictions are out of the scope of this 
work and, on the other hand, the reflected spectrum lies 
in the X-ray regime and is not expected to contribute 
much in the ISAF cooling (in which we are mainly inter- 
ested). Second, we kept the fraction of the flux that the 
disk reflects in the different phases of the oscillation con- 
stant. This, however, may not be the case (see for example 
Nayakshin et al. 2000). 

Keeping this in mind, we find that in the irradiation 
model q cxp — q CO nti < for the entire range of the parame- 
ter space r, T c , / that we have investigated. It seems to be 
a robust result that when soft photons are product of re- 
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processing of X-rays by the surface of the thin disk, small 
amplitude radial oscillations of the ISAF are damped. 

5.2. Proton illumination 

In the proton illumination model the soft photon flux 
comes from the overlap region only. Part of the proton illu- 
mination energy (about half, like in the Haardt-Maraschi 
energy balance) is reprocessed into blackbody radiation 
and the rest is emitted in the X-ray band due to the 'warm' 
layer that forms in the overlap region (see Deufel & Spruit 
2000 for details). Here we will ignore the X-rays coming 
from the warm layer. The total soft luminosity is 

L soft ~ ^F P A, (20) 

where F P is given by Eq. (JTJ and A — 2n(R 2 — R% T ) is 
the overlap surface area (the factor 2 accounts for the two 
sides of the disk) . Taking the Lagrangian derivative of Eq. 
(|20fl . we arrive at the result 

AL soft = 2i soft f _ jasm(u)t), (21) 

where K = Rtr/Rout, or / = 1 — K . From this equa- 
tion it follows that for K 2 < 2/3 (which corresponds to 
/ > 0.18), the soft flux decreases on expansion and in- 
creases on contraction of the ISAF. So for large overlap 
fraction, / > 18%, any small amplitude radial oscillations 
are expected to be damped (see also the discussion of the 
the synchrotron soft photon case in Sect. 3.1). Of more 
interest is the regime of / = 1 — Rt r /R on t < 10%. 

From calculations with different fractional overlap /, 
we find that there is a characteristic value f below which 
oscillations are excited. This is illustrated in Fig. 2. It 
shows, as a measure of the work integral, the quantity Y — 
(<7exp/<3coiitr — the fractional variation of the cooling 
rate, divided by the relative amplitude a of the oscillation. 
It is proportional to the growth rate r\ of the oscillations 
(see appendix A). For the example shown (r = 1.5 and 
T Q = 130 keV) the regime of excitation is / < 0.07. The 
results remain qualitatively the same for different choices 
of T and r. 

Until now we have only dealt with small amplitude 
oscillations. What limits the growth of these oscillations? 
To answer this questions we have also calculated the vari- 
ation of the cooling rate for larger amplitudes. The value 
of the amplitude at which the work integral changes sign 
is a good estimate of the amplitude reached by the oscilla- 
tion. Again we approximate the work integral by the sum 
of the integrand in the expanded and contracted phases. 

Fig- 01 (upper plot) shows the results for an equilib- 
rium overlap fraction / = 0.04. The growth rate has a 
maximum for a — 0.04. When a > 0.04 the ISAF, during 
part of the contracting phase, does not interact with the 
disk (/ = 0) and X-ray flux ceases. Note that the work 
integral is positive for all amplitudes: there is no limit 
to the growth in this simplified case. More realistically, 
hydrodynamic nonlincaritics ignored here would limit the 
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Fig. 2. The growth/damping rates of small amplitude 
oscillations for different overlap fractions / = (i? ut — 
i?tr)/-Rout in the proton illumination case. The left verti- 
cal axis shows the work integral Y = (q cxp — ?contr)/(<jt<Zeq)j 
the right axis the corresponding growth rate rj in units of 
the frequency uj of the oscillation. For / < 0.07 the work 
integral is positive and the oscillations are excited. 
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Fig. 3. Work integral Y of large amplitude oscilla- 
tions. Upper panel: the proton illumination case (as- 
terisks) and the case where there is also a soft syn- 
chrotron flux with Lm um /L S y n — 230 (diamonds). Lower 
panel: Lm um /L syn = 230 (diamonds), 163 (triangles), 150 
(squares), 138 (X). Note that for Li\\ um /L syn = 150, small 
amplitude oscillations are damped, while larger amplitude 
oscillations (a ~ 4%) are excited. 



oscillation. Also, other sources of soft photons will have a 
limiting effect, as discussed below. 

5.3. Combination of different soft photon sources 

The situation described in the previous section leaves out 
some physics. Synchrotron photons from the ISAF and 
soft photons produced by irradiation of the cool disk are 
expected to also play a role in determining the ampli- 
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Fig. 4. Work integral for the case where both proton il- 
lumination and irradiation are taken into account. Upper 
panel: small amplitude oscillations as a function of over- 
lap fraction /. For / < 0.05 small amplitude oscillations 
are excited. Lower panel: work integral as a function of 
oscillation amplitude, for fixed overlap fraction / = 0.04. 
In this case the oscillation grows up to a ~ 0.09. These 
results are for ISAF temperature ksT c = 150 keV, and 

T = 1. 



tude of the oscillation, since we have already shown that 
they act as damping mechanisms. As a first step in this 
direction, wc now include synchrotron photons into the 
simulation. The amount of synchrotron soft photons pro- 
duced by the ISAF depends on the unknown field strength. 
Hence we treat the synchrotron luminosity L syn as a free 
parameter. The soft flux due to proton illumination in this 
example is a black body with Tin um = 100 eV, while the 
typical energy of the synchrotron photons is set at ~ 3 
eV. The values of the other parameters are T c — 130 keV, 
r = 1.5, / = 0.04. 

For iiiium/isyn = 230 (Fig. 0, the oscillation ampli- 
tude grows up to a ~ 0.05 and stabilizes at that level. A 
weak synchrotron flux can limit the amplitude's growth to 
values slightly larger than the overlap fraction. However, 
if the synchrotron flux is strong enough, £in um /£ syn <140, 
the oscillation does not get excited at all. 

An interesting feature of the combined synchrotron- 
proton illumination model for the soft photons is the ap- 
pearance, in this example, of 'hard' oscillations. As Fig. [3 
shows, small amplitude oscillations are damped, while the 
work integral is positive for oscillations with larger am- 
plitudes, a ~ 0.04, if the synchrotron flux is in the range 
140 < L alum /L syn < 155. 

In this work we have ignored any variation of the vis- 
cous heating during the oscillation. However, viscous heat- 
ing can trigger/damp oscillations on the viscous time scale 
(e.g. Kley et al. 1993). Following the a viscosity prescrip- 
tion (Shakura & Sunyaev 1973) and taking that H ~ R for 
an ISAF, we have t v - mc ~ idyn/a;, where H is the thickness 
of the ISAF and idyn = l/u> is the dynamical time scale. 
So the viscous excitation/damping rate is r] v i sc /uj ~ a. 
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Comparing this with the rates due to radiative cooling 
(Fig. 2, 3), we see that for a < 0.1 the effects of viscous 
heating on excitation can be safely neglected. 

5.3.1. Reprocessed hard X-rays 

We now consider the case when the soft photon flux is 
produced by the joint effects of proton illumination and 
reprocessing of the hard X-rays from the ISAF. Keeping 
the proton illumination component of the soft flux as a 
black body with Tm um — 100 eV and, calculating the 
temperature profile of the disk due to X-ray irradiation 
with the method described in Sect. 5.1, we study whether 
small amplitude oscillations are excited or damped for dif- 
ferent overlap fractions /. In this example we have fixed 
k B T c = 150 keV and r = 1. 

The results are presented in Fig.|3| showing the growth 
rate as a function of the overlap fraction and oscillation 
amplitude. For / < 0.05 small amplitude oscillations are 
excited. Fixing the equilibrium overlap fraction / to 0.04, 
we have also calculated the variation of the cooling rate of 
the ISAF for larger oscillation amplitudes (lower panel). 
The maximum growth rate occurs at an amplitude a = 
0.04. The amplitude of the oscillation can be estimated 
from the change of sign of the grwoth rate as a function 
of amplitude, yielding a max ~ 0.09. 

5.3.2. Viscous dissipation in the cold disk as a source 
of soft photons 

A source of soft photons ignored so far is viscous dissi- 
pation in the thin disk. Including this soft-photon source, 
using standard thin disk expressions (Shakura & Sunyaev 
1973), we find that the results presented in the previous 
sections remain essentially the same. This is so because 
the intrinsic luminosity of the disk has only weak trigger- 
ing/damping effects on the oscillations of the ISAF. 

5.4. Spectral evolution during the oscillation 

From the discussion in the previous sections, it is evident 
that if the overlap fraction between the thin disk and the 
ISAF is small enough (typically smaller than 7%), radial 
oscillations of the ISAF can be excited. The overlap re- 
gion is a source of hard X-rays (Deufel et al. 2001, 2002). 
The ISAF itself, however, may actually be the main X-ray 
producing region, and its oscillations are expected to have 
direct implications on the shape of the emitted spectrum. 

On expansion of the ISAF the optical depth, the elec- 
tron temperature and hence the Compton y parameter 
drop. The soft photon input on the other hand increases. 
As a result the spectrum becomes softer while the total 
X-ray flux increases. 

The opposite happens in the contracting phase when 
the source becomes harder and the X-ray flux decreases. 
So the X-ray spectrum exhibits a characteristic pivoting 
behavior which is shown in Fig.0 In this figure (which cor- 
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Fig. 5. Spectra predicted by the combined model of Fig. 
3 for different phases of the oscillation. Solid line corre- 
sponds to the equilibrium position while dash-dotted lines 
correspond to maximum expansion and maximum con- 
traction. Ln\ nm / L syn = 230 (see diamonds in Fig. 3). 

responds to the case where both proton illumination and 
synchrotron photons are taken into account) we have plot- 
ted the spectrum computed by the Monte Carlo code in 
equilibrium position, in maximum expansion and in max- 
imum contraction. 

One clear conclusion from Fig. [5] is that the oscilla- 
tions are predicted to be soft i.e. to have more power in 
soft X-rays than in hard X-rays. A second conclusion is 
that a simultaneous oscillation (but of opposite phase) 
should be observed at optical/infrared wavelengths, if the 
synchrotron soft flux accounts for a large fraction of the to- 
tal optical/infrared luminosity of the source. Implications 
from these conclusions are discussed in the next sections. 

6. Consequences/Predictions 

Radial oscillations of the ISAF excited by the interaction 
between an ISAF and a cool disk have a number of char- 
acteristics from which they can be identified and conclu- 
sions can be drawn about the accretion geometry. The 
most promising source of excitation turned out to be the 
interaction between ISAF and cool disk via proton illumi- 
nation. For this mode of interaction, the model predicts 
the following. 

i) The QPO frequency is predicted to be close to the 
Kepler frequency in the outer edge of the ISAF. For fre- 
quencies characteristic of the oscillations (around 1 Hz), 
the ISAF must extend to i? ut ~ 10Q-R S . This is compa- 
rable to the extent of the ISAF in the low/hard state in 
Cygnus X-l, as deduced by Gilfanov et al. (2000) from 
Fourier frequency-resolved spectroscopy. 

ii) As shown in Sects. (5.2), (5.3) the excitation of the 
oscillations can be very fast, with tj/lu ~ 1. This implies 
that they are expected to show non-sinusoidal behavior 
and higher harmonics. Under these conditions our assump- 
tions about almost sinusoidal oscillations arc not valid any 
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more. The main conclusions, however, are expected to stay 
qualitatively the same. 

iii) There are cases with 'hard' or finite amplitude in- 
stability, in which small amplitudes are damped but finite 
amplitudes can grow rapidly. 

iv) If synchrotron flux accounts for part of the soft flux 
that is Comptonized, then infrared and/or visible mod- 
ulations of the same frequency with the QPOs are pre- 
dicted. The modulations will more likely be observed if 
synchrotron flux accounts for a large fraction of the total 
source emission in the relevant wavelengths. 

v) Excitation of the oscillations occurs for overlap frac- 
tion of the ISAF and the cool disk less than about 7%. If 
the observed QPOs are excited by the processes described 
here, there must be a rather sharp transition from the 
outer thin disk to the inner hot flow. 

A final conclusion is that the QPOs are predicted to 
be soft, showing larger rms amplitudes in the soft X-ray 
bands than in the hard ones. This, at first sight, seems to 
conflict with the usual observational trend that QPOs are 
hard. There are, however, several observations of sources 
in the low/hard state where soft QPOs have been ob- 
served. During the decay after the 2000 outburst of XTE 
J1550-564 (Kalemci et al. 2001), or in the hard state of 
Cygnus X-l (Rutledge et al. 1999) soft oscillations have 
been observed with frequencies around 1 Hz. This sug- 
gests that there are several mechanisms contributing to 
QPO excitation, but makes clear that other mechanisms 
than those described here need to be explored. 

7. Summary and discussion 

The implications of any model that attempts to describe 
QPOs depend strongly on the assumed accretion geom- 
etry. In this work we have assumed that the accretion 
flow consists of a standard thin disk which is truncated in 
a radius R tr and a hot two-temperature ion supported 
flow (ISAF) inside this radius. Due to its tendency to 
spread outward viscously, the ISAF extends to a radius 
-Rout > -Rtr, resulting in partial overlap of the two flows. 

We have focused on the frequency and the excita- 
tion/damping mechanisms of the fundamental mode of 
radial oscillations of the ISAF. The oscillatory frequency 
of this mode is close to the Kepler frequency in the outer 
edge of the ISAF. We find that the main contributor to the 
excitation of these oscillations is the variation of cooling 
rate of the ISAF in the different phases of the oscillation. 
When the ISAF cools more efficiently in the expanded 
phase than the contracted phase, an oscillation is excited. 

We have assumed Compton cooling to be the dominant 
cooling physical process of the ISAF. Since the source of 
the Comptonized 'seed' photons depends on the still un- 
known details of the accretion geometry, we have consid- 
ered three different physical mechanisms that have previ- 
ously been proposed to account for the soft photon flux: 
i) synchrotron flux due to frozen in magnetic fields in the 
ISAF, ii) irradiation of the thin disk by X-rays from the 



ISAF (Haardt & Maraschi 1993) and iii) proton illumi- 
nation of the disk in the overlap region (Deufel & Spruit 
2000). 

We have developed a Monte Carlo code to calculate 
the Compton scattering and to calculate the cooling of the 
ISAF, as well as the emergent X-ray spectra during the 
oscillations. We find that soft photons produced by syn- 
chrotron and by reprocessing of the hard flux act as damp- 
ing mechanisms of the oscillations. On the other hand, 
proton illumination excites oscillations if the overlap frac- 
tion / between the disk and the ISAF is small enough 
(of the order of 7%). Also models that combine proton 
illumination with synchrotron or irradiation soft flux can 
produce (if some constraints are fulfilled) radial oscilla- 
tions with amplitudes of the order of the overlap fraction, 
i.e. of the order 10%. 

The oscillations arc predicted to produce a QPO sig- 
nature in the X-ray power spectrum that has more power 
in the soft X-rays than in the hard X-rays (soft QPOs). 
This leads us to identify the low frequency QPOs ob- 
served in the decay of the 2000 outburst of XTE J 1550-564 
(Kalemci et al. 2001) and in the hard state of Cygnus X-l 
(Rutledge et al. 1999) within the framework of our model. 
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Appendix A: The decay/excitation time scale. 

The damping/growth rate rj of the amplitude of an oscil- 
lation can be denned as 



77 = 1/2 



(W) | 

{Eosc) 



(A.l) 



where () expresses the average value of a quantity over a 
period, E osc is the oscillatory energy and W is the rate of 
energy gain of the oscillation. If (W) is a negative number, 
then rj < and corresponds to the damping rate of the 
oscillatory amplitude, while if (W) is positive, then r\ > 
and corresponds to the excitation rate. 



The total energy of the oscillation equals the kinetic 
energy of the oscillation when the medium crosses the 
equilibrium position. With Eq. integrating over the 
volume of the ISAF we find 



5<7T 



(A.2) 



For nearly sinusoidal oscillations, (W) is given by the 
integral 



(W) = - 



T 



dm I (7 - 1) — Aqdt, 
P 



(A.3) 



where Aq is the variation of the cooling rate per unit mass 
of the ISAF. In general Aq is function of the space coordi- 
nates and time. Taking the cooling rate constant over the 
volume of the medium (numerical checks show that this 
is rather realistic) and to vary sinusoidally with time, we 
have 



Aq = q amp sin(ujt + <j>). 



(A.4) 



With <j amp = (g cxp -<7contr)/ (2Misaf) being the amplitude 
of the oscillation of q. Also <? cxp / <7 con tr is the cooling 
rate of the ISAF in maximum expansion/contraction and 
Misaf is the mass of the ISAF. 

Expression (|A.4(I contains a phase <fi. For — 7r/2 < <f> < 
7r/2, the oscillation is excited, while for tt/2 < <\> < 3/2tt, 
the oscillation is damped. We will focus on the extreme 
cases = and <j) — 7r. (The first corresponds to <? a mp > 
and the second to q a mp < 0). 

Integrating l|A.4|) over a period and over the mass of 
the ISAF we find 



(W) = 3(7 - l)a 



yexp ycontr 



(A.5) 



So, from (|A.lfl we can estimate the damping/growth rate 
of the amplitude of the oscillation. Of more interest is 
the quantity r]/u>, where to is the the frequency of the 
oscillation, taken to be the Keplerian frequency at the 
outer edge of the ISAF. 
This yields 



1 _ 15(7-1) 



L x (q, 



cxp ~ <7contr)/<Zeq 



4 t (GM/R out y/ 2 L Edd 



(A.6) 



where L x is the X-ray luminosity of the source and q eq 
is the the cooling rate of the ISAF in its equilibrium 
position. In order to derive this equation we have also 
used the expression for the Eddington luminosity Le = 
AnGMrripc/ <Jt and the fact that L x ~ g e q (since inverse 
Compton is responsible for both the cooling of the ISAF 
and the X-ray luminosity of the source) . 

The last expression can also be written more conve- 
niently as 



V 

id 



3-5 1/2 L x ((fcxp — <?contr) 



out r 



aq C q 



(A.l) 



Note that the growth rate of the oscillation is proportional 

to Y = (<j cxp - g CO ntr)/(a<7eq)- 



